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Abstract, y Doradus (y Dor) are F-type stars pulsating with high order g-modes. Their instability strip (IS) overlaps the red 
edge of the 8 Scuti (8 Set) one. This observation has led to search for objects in this region of the HR diagram showing 
p and g-modes simultaneously. The existence of such hybrid pulsators has recently been confirmed [ 10] and the number of 
candidates is increasing (e.g. [ 17]). From a theoretical point of view, non-adiabatic computations including a time-dependent 
treatment of convection (TDC) predict the existence of y Dor/5 Set hybrid pulsators (U, OH]). Our aim is to confront the 
properties of the observed hybrid candidates with the theoretical predictions from non-adiabatic computations of non-radial 
pulsations including the convection-pulsation interaction. 
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7 DOR/5 SCT HYBRID CANDIDATES 

There are presently three y Dor/5 Scuti hybrid pul- 
sator candidates, HD 49434 JH, HD 114839 01 and 
BD+18 4914 il. One more object, HD 8801 was al- 
ready proposed as y Dor/5 Scuti pulsator by 1112] and has 
recently been confirmed as a hybrid pulsator iflQfl . The 
available stellar parameters for these four stars have been 
collected from literature and summarized in Table 1 . In 
Figure[T]we plot their location in the HR diagram as well 
as the observational y Dor instability strip J 1111 and the 
red edge of the 8 Set instability domain 01911 . We note 
that these four stars have quite close T e ff (within 100 K), 
and are located near the blue edge of the y Dor IS and 
inside the 8 Set IS. 



THEORY VERSUS OBSERVATIONS 

To study the pulsation properties of these stars we have 
at our disposal a grid of stellar models computed with 
the evolution code CLES [22]. The grid properties are 
the following: stellar masses range from 1.2 to 2.5 M G 
with a step of 0. 1 ; four different chemical compositions 
described by the metal mass fraction Z = 0.01 and 0.02 
with a hydrogen mass fraction X = 0.70 and 0.73 are 
available. Moreover three different values of the mixing 
length parameter of convection (C£mlt = 1-4, 1.7, 2.0) 
and two values of the overshooting parameter (Oqv = 0.0, 
0.2) can be chosen. The pulsation analysis is done by us- 
ing a version of the non-adiabatic pulsation code MAD 



TABLE 1. Stellar parameters of the four hybrid y Dor/5 Set 
candidates. 



HD 8801 



HD 49434 



Spectral type 
Parallax n (mas) 
Teff (K) 
logg (cgs) 

[Fe/H] (dex) 
R(R e ) 

v sin i (km s ) 



A7m 
17.91 ±0.75 (lg) 



F1V 
24.94 ± 0.75 



7345 ± 155 



4.2 



(10) 



7300 ± 200 



(IS) 



(12) 
0.77 ± 0.03 



4.4 ± 0.2 



(23) 



(23) 



(10) 



0.825 ± 0.022 



(23) 



1.7 ±0.1 
55 ±5, 



(12) 
12) 



0.10 ±0.12, 
1.60 ±0.05 



87 ±4 



23) 
(16) 



(23) 



M (M s ) 


1.54 ± 0.03 (10 ) 


1.55 ±0.14 (2 ) 




HD 114839 


BD+18 4914 


Spectral type 


Am 


F5 (Am? (16) ) 


Parallax n (mas) 


5.04 ± 1.04(18) 




Teff (K) 


7356 ± 77 (16) 


7250 (19) 


logg (cgs) 


4.39 ±0.5 m 


3.77(i9) 


"*(£) 


1.132 ±0.18(5) 


0.92 (19) 


[Fe/H] (dex) 


0.04±0.15 (1) 




R(R e ) 


2.177 ±0.450 (16) 




v sin ( (km s ) 


66.7 ± 5.0 (1) 




M(Mq) 





that includes the effects of the convection-pulsation in- 
teraction (|0], J3]). In fact, it is necessary to include the 
effect of convection in order to match the observational 
red edge of the 8 Set IS and therefore to study the hy- 
brid pulsators. Dupret et al. Jit] found that a value of 
OCmlt = 2.0 is necessary to fit the location of observa- 
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FIGURE 1. HR diagram location of the four hybrid pul- 
sator candidates (from bottom to top: HD 8801, HD 49434, 
BD+18 4914, HD 114839) and their 1 a error boxes, except 
for BD+18 4914 whose T e g and L uncertainties were not avail- 
able. Dashed lines represent the observational y Dor IS 1 1 1 ] , the 
dotted-dashed line represents the red edge of the 8 Set IS fH^l . 
The empty square locates the model we selected for our non- 
adiabatic analysis. Evolutionary tracks computed with X = 0.7, 
Z = 0.02, GImlt = 2.0 and Ofey = 0.0 are shown in dotted lines. 

TABLE 2. Parameters of the H model 



M (M s ) 


1.54 


Kmlt 


2.0 


R{R Q ) 


1.552 


OJov 


0.0 


TeffOQ 


7346 


X 


0.70 


logg (cgs) 


4.24 


z 


0.02 




0.799 







tional 7 Dor and 8 Set IS (see also lfl3ll ). In this prelim- 
inary analysis, we have restricted our choice of parame- 
ters to X = 0.70, Z = 0.02, a ML T = 2.0 and a ov = 0.0. 
The theoretical instability domain is shown in Figure [2] 

We computed an additional main-sequence model 
(thereafter H model) whose HR diagram location is close 
to that of HD 49434 and HD 8801 (square symbol in Fig- 
ure[TJ. The parameters of that model are given in Table|2] 

The pulsation analysis of the H model for mode de- 
grees I = 0-8 reveals that it behaves as a hybrid pulsator, 
with unstable y Dor modes (from v = 0.789 c/d to 7.03 
c/d) as well as unstable 8 Set modes (from v = 1 1.0 c/d 
to 49.0 c/d) separated by a region of stable modes. To 
illustrate that we chose three different modes: a typical 
7 Dor g-mode (£=l, g24), a typical 8 Set /?-mode (£=l, 
P2) and a stable mode between / Dor and 8 Set frequency 
ranges (£=l, gi). In Figure [3] we present their propaga- 
tion diagram and their work integrals. Regions where the 
work increases (vs. decreases) are driving (vs. damping) 
the oscillations. For the y Dor g-mode, we see a clear 
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FIGURE 2. logr e ff-log.P theoretical instability domain for 
I = 0^1 modes for y Dor and 8 Set models with X = 0.7, 
Z = 0.02, a M LT = 2, Oov = and 1.2 M < M < 2.5 M Q . In 
grey are the y Dor and 5 Set domains and in black is the 
hybrid domain. 



driving mechanism at the location of the Fe opacity bump 
(logT ~ 5.3) but this K"-mechanism is not sufficient to 
globally excite the mode. The main driving occurs at the 
base of the convective envelope (CE) by the flux block- 
ing mechanism (FBM) (in agreement with J9J and (5|]). 
For the p-mode there is a contribution of the FBM at the 
base of the CE and a contribution from the fc-mechanism 
in the He partial ionization zone. For the stable mode, 
there is an efficient radiative damping in the inner layers 
of the star. Furthermore, the amplitude of the eigenfunc- 
tion in the outer layers is small due to the presence of a 
large evanescent region before reaching the base of the 
CE, which inhibits the FBM (see also {(&). 



Rotational splitting: application to HD 8801 
and HD 49434 

Even if theory predicts hybrid pulsators, we should 
wonder if the high frequency modes detected in HD 8801 
and HD 49434 are really 8 Set modes or prograde g- 
modes moved to higher frequencies due to rotational 
splitting. It is well known that for modes with a pulsa- 
tion frequency (PF) a comparable to, or lower than, the 
rotational frequency i2, the Coriolis force term plays a 
major role in the equation of motion and the perturba- 
tive approach is no longer valid. Since y Doradus stars 
show low-frequency g-modes, the effects of the Coriolis 
force cannot be neglected even if the star does not rotate 
fast. Dintrans & Rieutord [3] showed that the perturba- 
tive treatment of rotation is no longer valid for y Dor 



< 3 with v eq = v sin i = 55 km s 1 IU2I1 . In this case, 




split g-modes are not sufficient to explain the highest ob- 
served frequencies. Therefore the spectrum of HD 8801 
can most probably be attributed to hybrid pulsations. 



CONCLUSION 

Using non-adiabatic computations including TDC treat- 
ment for models with (Zmlt = 2.0, we predict the ex- 
citation of both 7 Dor and 8 Set modes separated by a 
region of stable modes in models located in the region 
of the HR diagram where hybrid candidates have been 
detected. Moreover, from a comparison between theoret- 
ical excited frequencies including the first order effect 
of rotation and observed frequencies of HD 49434 and 
HD 8801, we emphasize that it is necessary to consider 
the effect of rotation on PFs case by case in order to char- 
acterize these candidates as hybrid pulsators or as y Dor 
stars with g-modes split by rotation. 



FIGURE 3. Propagation diagram (top) & work integrals 
(bottom) for 1=1 §24 » gi an d P2 modes of the H model. Top 
panel: thick dark lines: propagation regions, thick grey lines: 
evanescent regions for g-modes. Bottom panel: vert, continu- 
ous line: base of the CE; vert, dotted lines: Fe opacity bump 
(logT ~ 5.3) & He partial ionization zone (logT ~ 4.9). 



with rotation period smaller than « 3 days. Moreover, 
one should also take into account the effect of rotation 
on the mode excitation I2H1 . Nevertheless, in a first ap- 
proximation, we estimate the rotational splitting by using 
the perturbative approach at first order: 



with (7 bs the PF in the observer frame, Ob the PF in the 
corotating frame, m the azimuthal order of the mode and 
/3 the Ledoux constant lla l. 

Uytterhoeven et al. 12311 identified some of the ob- 
served modes of HD 49434 as I =3-8 prograde modes, 
and estimated the value of the equatorial velocity to be 
236 km s _1 . Even if we adopt as equatorial veloc- 



I'eq 



ity v sin i = 87 km s , the 1 = 6 modes split by rotation 
can reach values of the order of the highest observed fre- 
quency (12 c/d). Therefore, the observed frequencies can 
be explained either by a combination of y Dor and 8 Set 
type modes, or by the splitting of high degree g-modes. 
Present observations do not allow us to confirm the hy- 
brid nature of HD 49434. 

Handler |l0[] performed a frequency analysis for 
HD 8801 using ground-based (GB) photometry. No 
mode identification is available but due to the limita- 
tions of GB photometry, we chose to restrict our study to 
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